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Abstract 

Very light pseudoscalars can arise from the symmetry-breaking sector in many extensions of the 
Standard Model. If their mass is below 200 MeV, they can be long-lived and have interesting phe- 
nomenology. We discuss the experimental constraints on several models with light pseudoscalars, 
including one in which the pseudoscalar is naturally fermiophobic. Taking into account the strin- 
gent bounds from rare K and B decays, we hnd allowed parameter space in each model that may 
be accessible in direct production experiments. In particular, we study the photoproduction of 
light pseudoscalars at Jefferson Lab and conclude that a beam dump experiment could explore 
some of the allowed parameter space of these models. 
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I. INTRODUCTION 



In many extensions of the standard model, the electroweak symmetry-breaking sector 
includes additional weak doublets or singlets. New CP-even, CP-odd and charged scalar 
states may be present in the physical spectrum. The masses of these particles are typically of 
the same order as the weak scale, and fine-tuning is required to make them much lighter. An 
exception occurs if the theory possesses an approximate global symmetry: a CP-odd scalar 
may become a massless goldstone boson in the limit that such a symmetry is exact, and a 
massive state that is naturally light in the case where the symmetry is only approximate. 
We will henceforth refer to such CP-odd states as light pseudoscalars. 

The most familiar example of a light pseudoscalar is the axion This pseudo- 

goldstone boson arises in a two-Higgs-doublet model with a global symmetry that allows 
independent phase rotations of the two Higgs fields. The axion arises as a consequence of 
spontaneous symmetry breaking and is exactly massless in the absence of gauge interactions. 
The axion acquires a small mass due to the QCD anomaly, which breaks this global symmetry 
at the quantum level. 

In other models, a global symmetry may be broken more significantly by a small param- 
eter that appears explicitly in the Lagrangian. For example, consider the Higgs potential 
for two Higgs doublets |3(, with a $2 ^ — ^2 symmetry: 

V = ^f^^ 1 + A /l^^ 2 +A 1 (^^ 1 ) 2 +A 2 (^^ 2 ) 2 +A3^^i^^2+A 4 |^t^2| 2 +^-((^t^2) 2 +li-c.) 

(1.1) 

In the limit A5 — > 0, this potential has a (J(l) x U(l) symmetry in which each doublet 
rotates by an independent phase. The spontaneous breaking of the diagonal U(l) symmetry 
yields a goldstone boson that is "eaten" when the theory is gauged; the remaining U(l), 
which rotates each doublet by an opposite phase, yields a physical goldstone boson state. 
When A5 is nonvanishing, this pseudoscalar develops a mass given by m 2 A = —\§v 2 , where 
v = 246 GeV is the electroweak scale. In this paper, we will consider pseudoscalars with 
masses in the 100 — 200 MeV range, for phenomenological reasons explained below. This 
can be achieved by setting A5 equal to a small number that is comparable to a light fermion 
Yukawa coupling — a light pseudoscalar would then be no more or less unnatural than a 
muon or light quark. 
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Of course, one can construct models in which the light fermion Yukawa couplings arise 
only via higher- dimensions operators in a more complete high-energy theory. The Yukawa 
couplings are identified with powers of the ratio of a symmetry breaking scale to the cut off 
of the theory, and therefore can be naturally small. By analogy, the U(l) symmetry present 
in the A5 = limit of Eq. (|1.1J1 may be broken by a field r\ that acquires a vacuum expecta- 
tion value at some high scale and contributes to the term of interest only through Planck- 
suppressed operators. Given this dynamical assumption, one predicts that the pseudoscalar 
mass is of the order ({jf) /M*) n l 2 v , where n is a positive integer, and M* = 2 x 10 18 GeV is 
the reduced Planck mass. Interestingly, for n = 2, and (rj) ~ 10 15 GeV (the nonsupersym- 
metric GUT scale), one obtains a pseudoscalar mass of approximately 100 MeV. One can 
imagine a variety of high energy theories in which similar results are obtained. 

Our interest in pseudoscalar masses between 100 and 200 MeV is motivated by the pseu- 
doscalar decay length and production cross section. We hope to have both in optimal ranges 
for detection of the pseudoscalar in possible photoproduction experiments at Jefferson Lab. 
As far as production is concerned, existing direct searches yield bounds on the pseudoscalar 
couplings that are weakest in this mass range, and a wide variety of experiments 0, U, 0, |f| 
severely constrain the pseudoscalar couplings for masses below 100 MeV. On the other hand, 
if the pseudoscalars are produced in significant but not overwhelming numbers, we hope for 
a decay length that is long enough to clearly separate the pseudoscalar decay signal from 
possible mesonic backgrounds. Pseudoscalars with masses above 200 MeV decay rapidly into 
muon pairs with a branching fraction near 100%, making detection via a separated vertex 
impossible. Thus, the 100 — 200 MeV mass window seems particularly promising for the 
experimental search that we propose in Section IV. 

To proceed with our phenomenological analysis, we must decide on the pseudoscalar 's 
couplings to standard model fermions; the pattern of these couplings is in fact quite model- 
dependent. In the standard two-Higgs-doublet models, the pseudoscalar couplings are pro- 
portional to Yukawa matrices multiplied by a ratio of the vacuum expectation values v\ and 
v%. On the other hand, one can employ simple discrete symmetries to construct three-doublet 
models in which only two doublets couple to quarks and do not mix with a third doublet 
coupling to the leptons. In this case, the pseudoscalar in the quark-two-doublet sector is 
entirely leptophobic. An analogous three-doublet model with a lepton-two-doublet sector 
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yields a pseudoscalar that has no couplings to quarks and is, hence, hadrophobic. Such 
models illustrate the range of the possible, but are not particularly well motivated. A much 
more appealing possibility is that the pseudoscalar may have no direct couplings to quarks 
or leptons at all. Let us comment on the motivation for such a fermiophobic pseudoscalar 
in more detail. 

One could imagine a number of reasons why a pseudoscalar may have suppressed couplings 
to standard model fermions. The suppression could be parametric, as in the type-I two 
Higgs doublet model when tan (3 is taken large. On the other hand, the suppression could 
be geometric, as in extra-dimensional scenarios in which fields have wave functions that are 
localized at different points in an extra dimension. Let us focus on a concrete realization 
of this second idea. Consider an S 1 /Z 2 orbifold of radius R, with standard model matter 
fields located at the y = fixed point, and gauge fields in the 5D bulk. Here y is the extra- 
dimensional coordinate. Assume that there exists additional vector-like matter in complete 
SU(5) representations (to preserve gauge coupling unification) as well as a gauge-singlet 
scalar field S, all isolated at the y = nR fixed point. A spontaneously broken approximate 
global symmetry of the singlet potential leads to a light pseudoscalar state that couples 
directly to the exotic matter multiplets only. The geometry of this scenario prevents mixing 
between the ordinary and exotic matter fields, which communicate with each other only via 
gauge interactions in the bulk. The scale of compactification can be taken large enough so 
that the effects of Kaluza-Klein excitations are irrelevant to the low-energy theory. 

Given the simplicity of the fermiophobic singlet scenario described above, we will focus 
our discussion on light pseudoscalars in the two-Higgs-doublet models of type-I and II and in 
the fermiophobic singlet scenario. We comment on the other possibilities where appropriate. 
In Section II, we analyze the experimental constraints on the light pseudoscalar in the con- 
ventional two-Higgs-doublet models, placing particular emphasis on the bounds from K and 
B meson decays. In Section III, the fermiophobic singlet scenario is studied, and in Section 
IV we study the possibility of detecting pseudoscalars of either type in photoproduction 
experiments at Jefferson Lab. Section V contains our conclusions. 
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II. CONSTRAINTS IN TWO-DOUBLET MODELS 



As we have described in the previous section, light pseudoscalars can arise in two-Higgs- 
doublet extensions of the standard model. Two popular options exist in which a discrete 
symmetry is imposed to forbid tree- level flavor changing neutral currents [3]: In Model I, 
all of the fermions couple to a single Higgs doublet, but none to a second. In Model II, 
the charge Q = 2/3 quarks couple to one Higgs doublet while the Q = —1/3 quarks and 
the leptons couple to another. A third possibility is that all fermions couple to both Higgs 
doublets, without the restriction of any discrete symmetry. An ansatz is then employed 
to make tree-level flavor changing Higgs couplings sufficiently small 8]. However, in this 
case it has been shown that a very light pseudoscalar will still lead to unacceptably large 
flavor- changing neutral currents . 

The coupling of the pseudoscalar Higgs to fermions is of the form ——Xff^fA where 
v = 246 GeV and Xf = cot (3 for all fermions in Model I, and Xf = cot (3 (tan ft) for the 
Q = 2/3 quarks (Q = —1/3 quarks and leptons) in Model II. Here tan/5 is the ratio of 
vacuum expectation values of the two Higgs doublets, and is a free parameter. 

There have been numerous discussions of the bounds on a light pseudoscalar, most re- 
cently by Larios, Tavares-Velasco and Yuan In Model II, the combined bounds from 
the nonobservation of J/ty —>■ Aj and T — > A'-f force tan (3 to be close to 1, since the former 
decay implies tan f3 ^ 1 and the latter implies cot /? < 1 |3j ; theoretical uncertainties don't 
quite allow the model to be excluded. In Model I, both decays imply only that cot (3 SA. 
Bounds from i], rj and w decays also force tan /3 ~ 1 in Model II and cot (3 < 1 in Model I |llj . 
Bounds from g — 2 are in flux at the moment, but do not appreciably change these results. 
(In addition, the g — 2 bound is only valid if one makes a strong assumption that there are no 
other possible nonstandard contributions at one loop.) Bounds from b — > sx Ap, and A\, 
can all be avoided by constraining the neutral and charged scalar masses jlfj. Thus, we will 
consider two cases: Model II with tan/5 ~ 1 and Model I with cot f3 ^ 1. After reviewing 
the decay modes and decay lengths of the light pseudoscalar, we consider the bounds from 
K and B meson decays, which present the strongest constraints on these models. 



5 



A. Decay Modes 



For a pseudoscalar lighter than twice the muon mass, there are only two possible decay 
modes, A —>■ e + e~ and A —>■ 77. The decay width into an electron pair is given by 

r — -^^Hlf 2 - (21) 

For tan/3 = 1, this gives a decay length of 0.6 — 1.2 centimeters in the pseudoscalar rest 
frame, for M A ranging from 100 to 200 MeV. This result scales as tan 2 (3 in Model II and 
cot 2 {3 in Model I. 

The decay into two photons proceeds at one loop with the width 

\J2 f N c Q}X f \ 2 a 2 M A 



A— >77 



Q4lT 3 V 2 



(2.2) 



where N c is 3 for quarks and 1 for leptons and Qf is the fermion charge. This expression 
is valid if the mass of the fermion in the loop is much larger than the momentum in the 
deoa, When this is not the case then the e Xact expression gl ven in Ret, g Q should 
be used. Note that Eq. ()2.2j) is independent of the heavy fermion mass. For the top quark 
contribution alone, with tan {3 = 1, one obtains a decay length in the pseudoscalar rest frame 
of 30 centimeters for m A = 100 MeV. Note that if one considers all quarks and leptons except 
the first generation fields, then the decay width is increased by a factor of 16, which would 
correspond to a decay length of 2 centimeters. For tan/3 ~ 1, the branching ratio into 
photons is 10% for M A = 100 MeV and 40% for M A = 200 MeV. Thus, we see that typical 
decay lengths, for tan (3 = 1, are on the order of a centimeter. For Model I with small cot {3, 
this decay length is increased by a factor of tan 2 {3. These decay lengths will, of course, be 
increased by a relativistic factor if the pseudoscalar has a large momentum (as it does in 
.B-decays) . 



B. K decays 

It is has been long known that the strongest bounds on axion models come from the 
decay K ,A Q fi; one expects that the sam e process win significantly const r a,n the 
light pseudoscalar scenarios of interest to us here. While many early analyses (that did not 
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take into account the heaviness of the top quark) seemed to exclude the possibility of a 
light pseudoscalar in the standard two-doublet scenarios, more recent work suggests that an 
allowed window remains. It was pointed out by Grzadkowski and Pawelczyk that there are 
two contributions to the decay amplitude and that the sum may vanish for some choices 
of model parameters ^j]. The first is a direct decay contribution involving the top quark 
and charged Higgs bosons at one loop; the second is an indirect contribution following from 
mixing between the axion and the 7r°, 77 and the rf. We refer the reader to Ref. |l4j for the 
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FIG. 1: The branching ratio for Kl — > ir°A for two values of tan/3 as a function of the charged 
Higgs mass. We choose Ma = 150 MeV. The experimental bound is approximately 4 x 10~ 8 . 



full expressions. As an example, the amplitude for 
schematically as 



7r + A in Model I can be written 



X w cot (3 F(m K ,m n ,m A ,m ri ,m vl ) + cot (3 G((3,m top ,m H +,U C KM) 



(2.3) 



The first term depends only on meson masses and is due to the pseudoscalar mixing; X w is a 
chiral Lagrangian parameter that is fixed by the data to be |A„,| = 3.2 x 10~ 7 [3]. The sign 
of A can be determined by matching chiral Lagrangian amplitudes to electroweak results [3J 
and is negative (the imaginary part is proportional to the CP violating factor e and 
is thus negligible). The second term represents the direct, one-loop decay amplitude, and 
depends on the top mass, the charged Higgs mass, and on CKM angles. Specifically, the 
second term may be written 



K) 



(2.4) 
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where 

£ = -^^^ s ^>Jcot^(A 1 + cot 2 /3A 2 ) . (2.5) 
Here Ai and A 2 are functions of the top, charged Higgs and W masses and are given explicitly 



/iiaiiu 

f. y. 



in Ref. [l3j- Numerically, the first term of Eq. ()2.5|) is typically a few times 10 -11 GeV and the 
second is typically 1CT 9 GeV. However, Grzadkowski and Pawelczyk show that the second 
term changes sign as the charged Higgs mass varies from 50 GeV to 1000 GeV, and thus at 
some value the total amplitude vanishes. We have plotted their results for the Kl decay in 
Fig. 1, setting tan/? = 1 (so our results then apply to both Model I and Model II), and also 
tan (3 = 50 in Model I. Consideration of K + and K$ decays leads to qualitatively similar 
results. 
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TABLE I: The allowed ranges for the charged Higgs mass (in GeV) for Ks, Kl, K^, and B decays. 
The four ranges overlap for all tan (3 shown. 

/,From Fig. ^ we see that there is a very narrow region of parameter space in which the 
branching ratio is suppressed. We now must consider whether the experimental bounds on 
Ks, Kl and K ± decays can be satisfied simultaneously. The Higgs Hunters Guide refers 
to two experiments Q| that search for the decay chain K + — ► n + A, A — ► e + e~, and 
obtain upper limits on the ttA branching ratio of order 10~ 8 . However, it is important to 



point out that a region between rriA = 100 — 150 MeV remains unconstrained due to the large 
background from the standard decay K + — > 7r + 7r°, followed by 7r° Dalitz decays. Without 
precise vertex detection, this can not be distinguished from the pseudoscalar signal. In the 
particular case of Model I with large tan /3, the decay length increases by tan 2 (3, and can 
be several meters. The pseudoscalar would then escape the detector. In that event, bounds 
from K + — > 7r + nothing 20. 21|. which range from 10~ 7 to 10 -10 , would apply. Again, 



the weaker (9(10~ 7 ) bound applies to a mass interval between tua = 130 — 160 MeV, as 
a consequence of larger experimental backgrounds. On the other hand, the experimental 
bounds on the decay Kl — > 7i°A, are uniformly strong over the entire range of pseudoscalar 



masses 



221 ]. Fortunately, one can fine-tune the charged Higgs mass to avoid contradiction 



with both charged or neutral kaon decay bounds. In Table I, we show the required range of 
charged Higgs masses for K + , Kl and Ks decays. It has been assumed that the A mass is 
150 MeV, so that the tighter experimental bounds in charged K decays apply; if the mass 
is between 100 MeV and 150 MeV, these bounds are relaxed and the ranges for K + and 
Ks decays are much wider. For all values of tan j3 shown in Table [fl the allowed ranges for 
charged Higgs mass overlap and all the bounds can be satisfied with a single fine tuning. For 
Model II, in which tan/3 ~ 1, the charged Higgs mass must be tuned to approximately one 
percent precision, but in Model I with larger tan/3, relatively mild fine-tuning is sufficient. 
Thus, kaon decays cannot completely exclude the existence of a pseudoscalar in the 100 — 200 
MeV mass range. 



C. B decays 

In B decays into KA, the pseudoscalar will have a relativistic gamma factor of 12 — 24, 
depending on its rest mass. Thus, the decay length into electrons will be approximately 25 
centimeters (times tan 2 j3). Because of the larger CKM mixing with the top quark, the Higgs- 
top loop contribution to the amplitude generally dominates over the mixing term by a larger 
amount than in the case of kaons. A simple estimate illustrates that the branching fraction 
is potentially large: The loop term involves CKM factors that are comparable to those found 
in tree-level semileptonic decays, while the 167r 2 in the loop is partly compensated by the 
smaller two-body phase space. The resulting prediction has a shape very similar to that 
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for K decays in Fig. 1. Again, there is a narrow region of parameter-space where the rate 
vanishes, and this region matches the narrow region in K-decays. This is not surprising since 
the analog of Eq. (|2.5|) for B decays has the same functional dependence on the charged 
Higgs mass, up to an overall factor. One might hope that higher order corrections would 
separate the K and B decay allowed mass windows, but a one-percent effect would not be 
sufficient to alter our qualitative results. 

What are the experimental limits? Recently, the BELLE Collaboration published a 
value for the branching fraction for B — > Ke + e~ of 0.75 ± 0.2 x 10~ 6 (2^. Since this is in 
agreement with theory, a bound on new physics contributions of approximately 2 x 10~ 7 can 
be obtained. However, the BELLE analysis included a mass cut on the electron-positron 
pair of 140 MeV, to suppress background from photon conversions and 7r° Dalitz decays. 
Thus, the bound does not apply to the 100 — 140 MeV window. The CLEO Collaboration 
has searched for B ± — > K nothin g a nd B° — > K° s nothing decays, and obtains a bound 
on the branching ratios of 5 x 10~ 5 [24|. While this does cover the mass range in which the 
BELLE analysis does not apply, it is only relevant if all the pseudoscalars escape detection. 
For masses between 100 and 140 MeV, one can ask what fraction of the A's will escape the 
detector. For tan (3 = 50, the decay length will be over 10 meters and almost all of the A's 
would escape; the CLEO bound would then apply. In general, approximately e~ 4cot2 13 of the 
A's escape the detector, which is a barrel calorimeter of roughly a meter radius. The bound 
would then be weaker by this factor, or 5 x icr 5 e 4cot2/3 for the branching ratio. Using this 
experimental bound, we find the allowed charged Higgs mass range given in Table I. We see 
that the same fine-tuning needed (for tan (3 ~ 1) for kaon decays will automatically suppress 
the B-decay rate. 

We conclude that neither model I nor II can be definitively excluded from the bounds 
from B decay, although fine-tuning is needed if tan/3 ~ 1, as required in Model II. 



D. Leptophobic Pseudoscalars 

As noted in the introduction, it is simple to have a three Higgs model in which two of the 
Higgs doublets couple to quarks (with Model I or Model II couplings) and a third couples to 
leptons. If the third doublet does not mix with the others, the leptonic couplings of the light 
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pseudoscalar are eliminated. The K and B decays discussed in the previous two subsections 
will generally not be affected in such a model. However, the decay of the pseudoscalar will 
now be entirely into photon pairs and the lifetime will generally be 2-3 times larger than 
the usual case. Note that in the 100 — 140 GeV mass window, the stronger bounds from Kl 
decays and from CLEO will certainly apply (without any significant exponential correction 
for decays inside the detector). Again, these bounds can be evaded with a suitable fine 
tuning of the charged Higgs mass. 



III. FERMIOPHOBIC PSEUDOSCALARS 



We have seen in the previous sections that a pseudoscalar state in the 100 to 200 MeV 
mass range is consistent with the stringent bounds from K and B meson decays. However, 
in the conventional scenarios considered thus far, this result follows from an accidental zero 
in the decay amplitudes, as well as a willingness to accept fine tuning. In this section we 
consider another possibility, that the couplings of the pseudoscalar to matter are naturally 
suppressed. After discussing the experimental bounds, we argue that a natural place to 
search for such a state is in a low-energy photoproduction experiment, such as those possible 
at Jefferson Lab. We estimate the production rate and comment on the relevant discovery 
signal in Section IV. 

We have already stated the motivation for considering a pseudoscalar state that is light: 
it might be the would-be goldstone boson associated with a global symmetry that is only 
approximate. In the introduction, we outlined a plausible scenario with a singlet scalar and 
a vectorlike multiplet in a complete SU (5) representation, taken to be a 5 + 5 for simplicity. 

Since the exotic matter is vector-like, it can be made arbitrarily heavy and integrated out 
of the theory. This leads to nonrenormalizable interactions between the pseudoscalar and 
the standard model gauge fields. If Mp is the mass scale of the vector-like matter tp, and 
the pseudoscalar coupling is given by (iA\/ V2)^ 5 ^, then one obtains 

C = % e^ pa AF. v F pa (3.1) 

for the effective coupling of the pseudoscalar to two photons. Here q represents the electric 
charge of tp, and F^ v is the electromagnetic field strength. Note that this can be generalized 
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to any non-Abelian gauge group by replacing q 2 with the Casimir Tp (defined by Tr[T a T b ] = 
Tp5 ab ) and by summing over the field strengh tensors. For a fermiophobic pseudoscalar in 
the mass range of interest to us, the only possible decay is to two photons, and from Eq. (|3.1|) 
we obtain the decay width 

, 16 a 2 X 2 m 3 A 

If Mp is not far above the top quark mass, say 200 GeV, and A = 1, then one obtains a 
lifetime 



r(M f = 200 GeV) = 1.1 x 10~ 3 sec fM^T) . (3.3) 

For energies of a few GeV, typical of the photoproduction experiments that we will mention 
later, the pseudoscalar can travel a macroscopic distance before it decays. A pseudoscalar 
with a mass of 150 MeV and an energy of 3 GeV will have a decay distance of 160 centimeters. 

One might think that the scenario described above is relatively insensitive to the bounds 
from meson decays due to the weakness of the pseudoscalar's coupling to ordinary matter. 
However, the experimental bounds on the branching fraction of K or B mesons to 7r+ pseu- 
doscalar are so stringent that operators like Eq. (J3.1|) are potentially significant, even when 
they contribute only at one loop. Here we estimate the contribution to K — > tcA in order to 
constrain the parameter space of the model. We comment on the constraints from B decays 
at the end of this section. 

The operator with the largest potential effect on low-energy hadronic decays is the glu- 
onic version of Eq, (13. lj) . We use a chiral lagrangian approach to estimate the branching 
fraction of interest [l4j. First we represent the light pseudoscalar nonet via the nonlinear 
representation 

S = exp(2i7r// 7r ) (3.4) 
where f n = 93 MeV is the pion decay constant, and where 7r is the matrix of fields 



TV 



(3.5) 



\ EZ. i(K°-K°) i + i. / 

\ V2 2^ A L K s) ^3 + / 

Here we have ignored CP violation and expressed the neutral kaons in terms of their CP 
eigenstates. Also note that we have chosen to include the T]', so that S is an element of U(3) 
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rather than SU(3). The X field transforms simply under the chiral SU(3) symmetry 

£ -> UlEU R (3.6) 

leading to the usual lowest order effective Lagrangian 

f 2 1 
Co = ^Tr d^d^E + ^/VTr (MS f + SM f ) , (3.7) 

where M represents the light quark current mass matrix. However, Eq. (|3.7|) does not take 
into account the QCD anomaly, which relates the divergence of the axial current to the 
product of gluon field strength tensors G^G^. A possible method of incorporating this 
effect into the chiral lagrangian is to introduce the additional terms 

Canom = ^q(x) log + Cq(xf (3.8) 

where q(x) represents 

q{x) = -&* G ^ ■ (3 - 9) 
Under an axial U(l) rotation, the field £ is multiplied by an overall phase, and it is not hard 

to show that Co + C anom yields the appropriate divergence of the axial vector current [25] . 

Now, one may treat q(x) as an auxilliary "glueball" field, and remove it using its equation 

of motion. One then finds 

1 / % , det £ x 



25] 



£ ~ = -4^U l0g deIStJ • (3 - 10) 
This term determines the rf mass, and the parameter c can be chosen accordingly. If one 
now includes the pseudoscalar coupling to gluons, an additional term must be added to 
Eq. (|3.8jh namely Aq(x) / (2\^2Mp) , in which case Eq. f!3.10|) is modified 

1 / / . det S 1 



Canom = log r + = A . (3.11) 

4c V 2 S det St 2V2M F J K J 

This interaction leads to mass mixing between the pseudoscalar and the T]'; we find that the 
mixing angle is given approximately by 

^ " zjsm (3 - 12) 

or numerically, 7 x 10 -5 • (200 GeV/Mp). We may extract the AS = 1 Ktct]' vertex from 
the chiral Lagrangian term 

Cas=i = ^jTr (Xuhd^d^) (3.13) 
13 



where h is octet-dominant AS = 1 spurion 

/ 

h = 

V 




1 




(3.14) 



and \ w = 3.2 x 10 is a parameter that takes into account the strength of the weak 



interactions |l5|]. We find 



V{K + -> 7i + A) 



1 \ 2 fl 2 



3847T m\fl {mA 

As a point of reference, if one sets m A = 100 MeV, one obtains the branching fraction 
5.6 x 10~ 7 ■ (200 GeV/M F ) 2 . 

Different experimental bounds are relevant depending on the lifetime and boost of the 
pseudoscalar. If the pseudoscalar decays inside the experimental detector, the relevant bound 
on the K + branching fraction is [3J 



+ 2m 2 K f [{m 2 K -m 2 w + m 2 A f - Am 2 K m 2 A } 1/2 . (3.15) 



BF(K + -> 7T+77) < 1.4 x 10" 



(3.16) 



If the pseudoscalar escapes the detector unobserved, one must contend with more stringent 
bounds, ranging from ~ 10~ 7 to ~ 10" 10 , depending on the pseudoscalar mass [2l|. In Fig. [2] 
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FIG. 2: Allowed parameter space for the fermiophobic scenario. 

we display the allowed region of the model's parameter space. Within the two excluded 
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regions toward the top of the figure, the pseudoscalar is long lived enough to escape the 



detector, while the branching fraction exceeds the bounds given in Ref. |2l| . The gap 
between these regions corresponds to a mass interval in which there are larger experimental 
backgrounds. Immediately below each of these excluded regions, the pseudoscalar decays to 
two photons within the detector (assumed to have a fiducial length scale of 1.45 meters 26^ 
and the weaker bound in Eq. ()3.16|) becomes relevant. However, one never reaches the region 
of parameter space excluded by the K + — > 7r + 77 bound since the vector-like matter would 
itself become light enough to be detected in direct collider searches. We will restrict ourselves 
to the allowed regions of Fig.|2]with smallest Mp in discussing pseudoscalar production rates, 
in the next section. 

Finally, we should comment on the bounds from the analogous decays of neutral kaons 
and B mesons. First, the K® indeed may decay into ir° A; however, the total width of the 
is approximately two orders of magnitude larger that that of the K + , so the branching 
fraction to the decay mode of interest is suppressed by this factor relative to our previous 
results. We therefore obtain no further bounds. The K^, on the other hand, has a total 
width that is about a factor of four smaller than that of the charged kaon. However, the 
decay K\ — > tt°A is CP violating, so that the decay amplitude is suppressed by an additional 
CP- violating spurion factor of ~ 10~ 3 [16[, and again no further bound is obtained. In the 
B system, the decay B — > Krj' is observed, and has a branching fraction of order 10~ 5 [27]. 
Using our previous result for the Arj mixing angle, we estimate that the branching fraction 
for B — > KA is O(10~ 15 ) and no further bound is obtained. 



IV. PRODUCTION AT JEFFERSON LAB 



We have seen that there is a window for light pseudoscalars in the 100 — 200 MeV mass 
range. For the two-doublet Model II (or Model I with tan/? ~ 1) the window requires 
substantial fine-tuning of the charged Higgs mass; for the two-doublet Model I with large 
tan p, there is less fine-tuning, and for the fermiophobic case there is a very large region of 
allowed parameter space. How can one detect these pseudoscalars? 

A number of authors have considered light pseudoscalar detection at high-energy collid- 
ers [ll], 13]. Larios, Tavares-Velasco and Yuan [l^] discussed production at the Tevatron, 
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the LHC and future colliders. They focused on the two-photon decay mode, which at high 
energies registers as a single photon signature. In this section, we consider the possibility of 
detecting the pseudoscalars we have discussed in a beam dump experiment at the Thomas 
Jefferson National Accelerator Facility (Jefferson Lab). 

Jefferson Lab has a high intensity photon beam directed into the CLAS detector in Hall 
B. The maximum energy is currently 6 GeV with an upgrade to 12 GeV planned. The pho- 
ton beam has a bremstrahlung spectrum with a luminosity of approximately 10 34 cm -2 sec -1 
if the photons are untagged. At the 12 GeV upgrade a monochromatic 9 GeV photon beam 
will also be available, with a luminosity of approximately 10 33 cm -2 sec -1 . The amplitude 
for pseudoscalar photoproduction may receive two possible contributions. In the conven- 
tional two-Higgs-doublet models, one can photoproduce the pseudoscalar most copiously 
off the strange quark sea in the proton. Second, in all the models we have discussed, the 
pseudoscalar may bremstrahlung off the incident photon via the loop-induced A77 vertex. 
Once produced, the pseudoscalar will travel some distance and then decay into either e + e~ 
or 77, depending on the model. If the beam dump consists of a meter or more of material, 
then most of the 77 background events will be suppressed. It is thus important that the 
lifetime of the A be sufficiently long that a substantial number make it through the beam 
dump. 

We first concentrate on production. Consider photoproduction of the pseudoscalar off 
the strange quark in the proton. The parton level cross section in the center of mass frame 
is 

da h 2 e 2 p 



dcosO 1447rs 3 / 2 



m 2 — i 2m 2 m 2 , s — m 2 

j b rv j 



s - m 2 (s - m 2 s ) 2 m 2 - t 

(4.1) 



2m 2 s m 2 A ts — {m\ + m 2 )(s + t) + m\ + m\ 



(m 2 -t) 2 (s-m 2 )(m 2 -£) 
Here, h is the Yukawa coupling of the A to the strange quark, p is the A momentum; we 
have approximated the initial photon momentum as V§/2 in the phase space factors to 
simplify the expression. In finding the full cross section for photoproduction, we multiply 
by the parton distribution function for the strange quark and integrate. However, since the 
parton model becomes less reliable at small momentum transfers, one must keep in mind 
that there is significant theoretical uncertainty from the small x region of integration, where 
the partonic cross section is largest. We therefore cut off the x integration at a value where 
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s — xs — 1 GeV 2 . We believe that this choice is reasonable. At lower s there will not be 
enough energy to produce 0, r\ and K mesons, and thus one expects an additional suppression 
from the electromagnetic form factor due to the decrease in available exclusive channels. The 
resulting cross section is rather insensitive to the beam energy, varying from 3.6 cot 2 (3 to 
2.0 cot 2 j3 femtobarns as the photon energy varies from 2 to 12 GeV. For a luminosity of 
10 34 cm -2 sec -1 , this will yield approximately 800 cot 2 (3 events per year. In order to be 
detected, these pseudoscalars must travel through a beam dump. The lifetime, as discussed 
in Section II, gives a decay length for a 100 MeV pseudoscalar of 1.2 tan 2 (3 centimeters times 
the relativistic factor of E/M A . Consider a 6 GeV beam and tan (3 = 1. The decay length 
is then 72 centimeters, and roughly 25% of the particles, or 200 particles/year, will travel 
through a one-meter beam dump. Since the differential cross section has a t-channel pole in 
the massless quark limit, it is forward peaked and this estimate will not suffer a substantial 
solid angle dilution. As tan/3 increases, the production cross section drops, but the decay 
length increases. In Table II, we show the number of events that traverse a one-meter beam 
dump per year, assuming 10 34 cm -2 sec -1 luminosity. These pseudoscalars will primarily 
decay into an electron-positron pair. One should keep in mind that the uncertainties caused 
by the low x cutoff could be substantial, and thus these event rates are approximate. Also, 
for larger beam energies and larger tan (3, the decay length will be too long for a substantial 
number of events to occur in a detector. Nonetheless, the relatively high event rate indicates 
that further experimental analysis is warranted. 

The second production mechanism is through the A77 vertex. In the two-Higgs-doublet 
models, the production mechanism already considered strongly dominates, but in the fermio- 
phobic model, pseudoscalar bremstrahlung off the incident photon is the only possibility. The 
parton level cross section is 

da X 2 Q 2 e 6 p{2m 2 q m\ + t 3 - 1[m\ - s)i 2 + ((m\ - m 2 q f + m 4 + 2s 2 - 2(m\ + 2m 2 )s)i) 

JcosO ~ 20 487r 5 M|t 2 s 3 / 2 

(4.2) 

where A is the coupling of the fermion in the loop to the A, M F is the mass of the fermion 
in the loop, p is the final state 3-momentum of the A, and Q is the quark charge in units of 
e. In deriving this expression, we have assumed that Mp is much greater than the photon 
energy (certainly true for the fermiophobic case). For \f% >> m q + m^, we find that a is 
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tan/3 4 GeV 6 GeV 9 GeV 12 GeV 24 GeV 



1 


125 


210 


350 


340 


340 


2 


150 


147 


144 


135 


120 


3 


88 


85 


77 


66 


46 


4 


53 


53 


48 


40 


32 


5 


41 


34 


28 


25 


21 



TABLE II: The number of pseudoscalars traversing at least one meter for various values of the 
beam energy and tan (3 in the two-doublet model with photoproduction off the strange quark sea. 
We have assumed a luminosity of 10 34 cm -2 sec -1 and a pseudoscalar mass of 100 MeV. Most will 
decay into an electron-positron pair. 



well approximated by 



a 



a'Q'X 2 1 
64vr 2 M 2 § 



(2s -m 2 A y log 



js~rn 2 A ) 2 
m 2 m 2 A 



3(s-m 2 A ) 2 



(4.3) 



The exact parton-level total cross section is shown in Fig. El The approximate expression 



0.15 - 



t 0.1 -, 




0.05 



4 6 
COM energy (GeV) 



FIG. 3: Parton-level production photoproduction cross section in the fermiophobic scenario as a 
function of center of mass energy, with Mp = 200 GeV and A = 1. 



given in Eq. (J4.3|) yields results that are visually indistinguishable from those shown in 
Fig. El Using CTEQ set 5L structure functions for the up and down sea and valence quarks 
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we obtain the total production cross section shown in Fig. |3J Assuming a monochromatic 
photon beam and a Jlab-like luminosity of 10 34 cm -2 s _1 , one estimates 315 production 
events per year per femtobarn of total cross section; qualitatively speaking, Fig. |U suggests 
O(10 2 ) events per year at an energy-upgraded Jlab, or at some similar facility. 

0.55 

0.5 
0.45 

0.4 

& 0.35 

c 

0.3 
0.25 

0.2 
0.15 

FIG. 4: Photoproduction cross section in the fermiophobic scenario, as a function of photon beam 
energy in the lab frame, with Mp = 200 GeV and A = 1. 

A more realistic analysis would take into account that the highest luminosity photon beam 
at Jlab is not monoenergetic, but has a bremstrahlung spectrum. We approximate this effect 
by assuming a total luminosity of 10 34 cm -2 s _1 , with a distribution dLjdE 1 oc 1/E 7 , with 
Ej ranging from 1 GeV up to the beam energy; the event rate is determined by the integral 

/ °IE dE -< ' < 4 ' 4 > 

Table III shows the events per year for a number of different choices for the beam energy 
and pseudoscalar mass. Unlike the two-doublet model, the lifetime discussed in Section 
III is sufficiently long that most of these pseudoscalars will traverse a one-meter beam 
dump. Another major difference is that these pseudoscalars will decay into two photons, 
i.e. they will look like long-lived 7r°'s. A more detailed analysis taking into account possible 
experimental acceptances and cuts would be needed to determine whether this signal could 
be separated from background under realistic conditions. 
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E 1 (GeV) m A = 100 MeV m A = 200 MeV 



6 



99 



79 



12 



108 



87 



24 



117 



96 



TABLE III: Photoprodution event rate per year in the fermiophobic scenario, with Mp = 200 GeV 
and A = 1. The total luminosity is taken to be 10 34 cmr 2 s~ l and a Bremstrahlung photon 
spectrum is assumed between 1 GeV and the beam energy. 

V. CONCLUSIONS 

We have considered light, elementary pseudoscalars with masses between 100 and 
200 MeV. We have argued that such states may evade the stringent bounds from K and B 
meson decays, while remaining of interest in searches at low-energy photoproduction exper- 
iments, such as those possible at Jefferson Lab. In conventional two-Higgs doublet models, 
light pseudoscalars may evade the strange and bottom meson decay bounds due to a possible 
cancellation in the decay amplitude. In this case, the coupling of the pseudoscalar to quarks 
is substantial and one can produce the pseudoscalar state copiously via photoproduction off 
the strange quark sea in a nucleon target. On other hand, if one wishes to avoid fine tuning 
in evading the decay bounds, one can consider very natural scenarios in which the pseu- 
doscalar is fermiophobic. We have presented one concrete realization of this idea, motivated 
by extra dimensions, and have isolated the allowed parameter space of the model. In the 
fermiophobic scenario, the pseudoscalar-two photon coupling leads to production via pseu- 
doscalar bremstrahlung off the incoming photon line. The event rate is substantial enough 
to make accelerator searches of potential interest. 
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